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ABSTRACT 


In  this  paper,  the  problem  of  3-dimensional  guidance  law  design  is  considered  and  a  new  guidance  law 
based  on  partial  sliding  mode  technique  is  presented.  The  approach  is  based  on  the  classification  of  the 
state  variables  within  the  guidance  system  dynamics  with  respect  to  their  required  stabilization 
properties.  In  the  proposed  law  by  using  a  partial  sliding  mode  technique,  only  trajectories  of  a  part  of 
states  variables  are  forced  to  reach  the  partial  sliding  surfaces  and  slide  on  them.  The  resulting  guidance 
law  enables  the  missile  to  intercept  highly  maneuvering  targets  within  a  finite  interception  time. 
Effectiveness  of  the  proposed  guidance  law  is  demonstrated  through  analysis  and  simulations. 
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1.  Introduction 

Design  of  guidance  laws  for  intercepting  of  highly  maneuver¬ 
ing  targets  is  a  difficult  nonlinear  problem.  The  missile  must 
intercept  its  maneuvering  target  within  a  finite  interception  time, 
in  some  case,  only  in  a  few  seconds.  Recently,  nonlinear  control 
theories  have  been  used  in  design  of  robust  guidance  laws.  Laws 
such  as:  Lyapunov-based  nonlinear  guidance  laws  [1-3],  first- 
order  sliding  mode  guidance  laws  [4-8],  and  nonlinear 
guidance  laws  [9-11],  are  obtained  based  on  Lyapunov  theorems 
on  asymptotic  stability  or  exponential  stability  of  all  states.  In 
[12,13],  guidance  laws  are  obtained  by  applying  a  second-order 
sliding-mode  control  method.  However,  these  laws,  because  of 
complexity  in  the  structure,  only  have  been  applied  to  planner 
interception  geometry. 

It  was  shown  in  [14]  that  in  a  practical  situation,  asymptotic 
stability  behavior  is  not  a  realistic  behavior  for  all  states  of  the 
guidance  system  and  a  new  guidance  law  based  on  partial  stability 
principle  was  designed  for  planner  interception  geometry.  The 
approach  was  based  on  the  classification  of  the  state  variables  of 
the  guidance  system  into  two  parts  with  respect  to  their  required 
stabilization  properties  where  stability  is  desired  behavior  for  only 
the  first  part  of  state  variables. 

In  this  paper,  a  new  3-dimensional  missile  guidance  law  based 
on  Partial  Sliding  Mode  technique  (PSM)  is  designed.  The  PSM 
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technique  was  first  introduced  in  [15]  and  in  this  work  it  is 
applied  to  the  guidance  control  problem. 

Most  of  the  existing  methods  in  the  field  of  nonlinear  control 
theory  are  based  on  stabilization  of  all  states  and  in  spite  of 
applications  of  partial  stability  in  many  of  engineering  fields  [16], 
there  are  a  few  papers  regarding  the  design  of  control  laws  which 
stabilize  only  a  part  of  system’s  states.  Also,  among  the  existing 
papers  in  the  field  of  partial  control,  most  of  them  only  consider 
a  case  study  and  try  to  design  control  laws  for  partial  stability  of  their 
specific  applications  [17-19]  and  do  not  pose  a  general  framework  in 
designing  partial  control  laws.  It  is  worth  noting  that  the  control 
schemes  posed  in  Refs.  [17,20-22]  are  uneasy  to  realize  and  are 
usable  only  for  systems  with  some  special  structures.  The  PSM 
technique  leads  to  design  of  a  partial  stabilizing  controller  for  a 
nonlinear  system  in  the  presence  of  uncertainties  and  external 
disturbances.  In  this  method  partial  sliding  surfaces  are  designed 
such  that  restricting  the  motion  on  these  surfaces  guarantees 
asymptotic  stabilization  of  only  the  first  part  of  state  variables.  This 
is  in  contrast  to  standard  sliding  mode  (SM)  technique  that  stabilizes 
all  state  variables  by  forcing  them  to  reach  sliding  manifolds  (which 
are  functions  of  all  states)  and  slide  on  them  toward  the  equilibrium 
point.  Another  structural  difference  between  PSM  and  SM  is  that  in 
PSM,  the  reduced  input  vector;  was  designed  while  in  the  SM,  the 
input  vector  is  wholly  designed.  This  idea  makes  the  possibility  of 
converting  the  control  problem  into  a  simpler  one  having  fewer 
control  input  variables.  Additionally,  the  remaining  control  variables 
make  the  possibility  of  imposing  some  constraints  on  the  behavior  of 
the  other  state  variables  [15]. 

In  the  present  work,  the  PSM  method  was  applied  in  guidance  law 
design.  This  approach  is  advantageous  from  practical  view-points, 
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since  it  forces  the  states  of  the  guidance  system  to  behave  as  it  is 
desirable  in  practice  and  is  compatible  with  a  successful  guidance 
scenario.  In  addition,  the  acceleration  vector  of  the  target  is  assumed 
as  an  external  bounded  disturbance  and  only  its  bound  is  required  in 
the  design  of  the  guidance  law  and  the  accurate  measurement  of 
target  acceleration  during  the  flight  is  not  necessary.  Effectiveness  of 
the  designed  guidance  law  in  achieving  zero-miss-distance  within  a 
finite  interception  time  is  demonstrated  analytically  and  also  through 
computer  simulations. 

2.  Problem  definition 

2.1.  Missile/ target  kinematics  model 

In  this  section,  a  3-dimensional  air  to  air  engagement  is 
considered.  The  vehicles  are  modeled  as  point  masses.  The  relative 
motion  between  the  missile  and  the  target  is  described  in  the 
spherical  coordinate  system  (r,0,</>)  with  the  origin  fixed  on  the 
center  of  the  mass  of  the  missile.  The  3-dimensional  missile-target 
geometry  is  shown  in  Fig.  1,  where  r  is  the  relative  distance 
between  the  missile  and  the  target.  Also,  6  and  </>  are  yaw  and 
pitch  line  of  sight  angles,  respectively.  Let  (er,e0,e^)  be  unit  vector 
along  the  coordinate  axes  and  =  r~er  be  the  relative  distance 
vector  along  the  line  of  sight  (LOS)  vector  (i.e.,  the  vector  connect¬ 
ing  the  missile  to  the  target  at  each  moment  with  direction  toward 
the  target). 

The  differentiation  of  gives  the  3-dimensional  relative 
velocity  as 

~r  =  r~er  +  vb  coscj)~e  e  +  r(!)~e  ^  (1) 

where  r  is  the  radial  component  of  relative  speed.  Also,  r6  cos 0 
and  r0  are  tangential  components  of  relative  speed  between  the 
missile  and  the  target,  respectively.  Differentiating  Eq.  (1)  yields 
the  components  of  the  relative  acceleration  as 

{•2  -2 

r-r0  -rO  cos2  </>  =  ojr-ur 

r6  cos  4>  +  2 r6  cos  </>-2 r<j)0  sin  </>  =  coe-ue  (2) 

r0  +  2f(/>  +  r02  cos  </>  sin  0  =  co^-U0 

where  co  =  (cor,co0,co^)r  and  u=(ur,ue,Uf))T  are  the  acceleration 
vectors  of  target  and  missile,  respectively.  Here,  acceleration  of 
the  target  is  assumed  as  an  external  bounded  disturbance.  It  is 
worth  noting  that  only  this  bound  is  required  in  the  design  process 
of  guidance  law  and  the  accurate  measurement  of  target  accelera¬ 
tion  during  the  maneuver  is  not  necessary.  In  the  sequel,  we 
consider  x(t)  =  [0(t),(/>(t),0(t),(/)(t),r(t),f(t)]r  as  the  state  vector. 


Z 


Fig.  1.  3  -dimensional  missile-target  geometry. 


Note :  Usually,  the  initial  conditions  of  engagement  are  such 
that  r0  >  0  and  r0  <  0  (as  assumed  in  the  present  research).  It 
means  that  the  target  is  in  front  of  the  missile  and  the  missile  is 
approaching  it. 

2.2.  Desirable  behaviors  for  each  state  variable 

For  the  interception,  it  is  sufficient  that  r(t)  becomes  zero  in  an 
instance,  say  tf,  where  tf  is  the  interception  time,  i.e.,  fitf)= 0  (zero- 
miss-distance).  Clearly,  it  is  not  convenient  for  r(t)  to  asymptoti¬ 
cally  converge  to  zero.  In  other  words,  the  asymptotic  convergence 
is  not  an  ideal  behavior  for  r(t).  It  should  be  noted  that  asymptotic 
convergence  behavior  means  that  the  missile  approaches  the  target 
within  an  infinite  time  interval.  It  is  evident  that  such  a  behavior  is 
not  a  desirable  behavior  in  missile  guidance,  which  is  supposed  to 
hit  the  target  with  a  nonzero  speed  in  a  finite  time.  For  this 
purpose,  it  is  sufficient  that  the  relative  radial  speed  between  the 
missile  and  the  target  satisfies  the  following  proposition. 

Proposition  1.  In  order  to  intercept  the  target  within  a  finite  time 
interval  {fitf)  =  0),  it  is  sufficient  to  have : 

e  [t0  tf)  s.t.  r  <  -C  <  0  Vt  e  [t!  tf]  (3) 

where  f  is  a  positive  constant  Also,  r(t)  is  positive  and  continuous. 
Consequently,  to  reach  a  zero  value,  i.e.  zero-miss-distance,  within  a 
finite  time,  it  should  decrease  from  its  value  at  tfi  <  tf)  down  to  a  zero 
value  at  tf.  Inequality  (3)  implies  r(t)-r(ti) < -C(t-ti),  which 
shows  that  fit)  reduces  within  the  time  interval  [fi  t].  Since  it  is 
desirable  to  have  fitf)= 0,  thus  this  leads  to 

,4, 

Hence,  our  expectation  of  r  is  that  it  satisfies  Proposition  1. 
In  other  words,  convergence  and  stability  of  this  state  is  not  under 
consideration.  Although  regulating  r  to  a  negative  constant  c  guar¬ 
antees  interception  of  the  non-maneuvering  target  within  a  finite 
time,  it  is  not  efficient  for  interception  of  a  highly  maneuvering  target 
in  an  acceptable  interception  time.  In  such  a  case,  to  improve  the  per¬ 
formance,  c  should  be  time  varying;  however,  for  its  determination, 
the  maneuvering  model  of  the  target  should  be  known  for  the  missile, 
while  it  is  not  always  possible  in  practice  [14]. 

Regarding  the  other  state  variables,  i.e.,  6,  6,  </>,  </>,  in  contrast 
to  two  first  states,  i.e.,  r  and  r ,  asymptotic  convergence  behavior  is 
desirable.  An  appropriate  guidance  law,  in  addition  to  decreasing 
the  relative  distance,  must  keep  the  pitch  and  yaw  LOS-angular 
rates  as  small  as  possible  [1 1  ].  It  means  that  it  is  desirable  to  have 
(b,  </>)  ->(0,0),  0->Ci  and  </>-»c2,  where  Ci  and  c2  may  be  free  or 
some  pre-specified  constants.  Therefore,  the  state  vector  may  be 
separated  into  Xi  =  [6  <f  6  <f]T  and  x2  =  [r  r]T  where  the  asymptotic 
stability  behavior  is  desirable  only  for  Xi. 

3.  Partial  sliding  mode  control  method 

In  this  section,  a  robust  partial  stabilization  technique  for  a  class 
of  nonlinear  dynamical  systems  based  on  the  first  order  sliding 
mode  control  idea  is  presented  [15].  This  partial  control  method 
makes  the  possibility  of  converting  the  control  problem  into  a 
simpler  one  having  fewer  control  input  variables  while  in  other 
existing  papers  in  the  field  of  partial  control,  the  input  vector  is 
wholly  designed.  For  this  purpose,  the  state  vector  (xeRn)  is 
separated  into  two  parts:  XieRni  and  x2eRni,  (rq- \-n2=n )  and 
accordingly  the  nonlinear  dynamical  system  is  divided  into  two 
subsystems.  The  subsystems,  hereafter,  are  referred  as  Xi  -subsys¬ 
tem  and  x2 -subsystem.  The  reduced  control  input  vector  (the 
vector  that  includes  components  of  input  vector  appearing  in  the 
Xi -subsystem)  is  designed  based  on  the  new  sliding  mode  control 
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technique,  i.eM  “partial  sliding  mode ”  in  such  a  way  to  guarantee 
asymptotic  stability  of  the  nonlinear  system  with  respect  to  the  first 
part  of  states  vector,  i.e.,  Xi  for  every  x2.  Furthermore,  the  remaining 
control  variables  (which  do  not  appear  in  the  x,\ -subsystem  but 
appear  in  the  x2 -subsystem)  make  this  possibility  to  impose  some 
constraints  on  the  behavior  of  the  second  part  of  states  (x2).  In  brief, 
there  are  two  main  structural  differences  between  partial  sliding 
mode  and  the  first  order  sliding  mode  techniques.  First,  in  partial 
sliding  mode,  sliding  surfaces  is  a  function  of  only  Xi ;  while  in  the 
conventional  first  order  sliding  mode  technique,  the  sliding  surfaces 
is  a  function  of  x.  Second,  in  partial  sliding  mode  method,  the 
reduced  input  vector  is  designed  while  in  the  conventional  one,  the 
input  vector  is  wholly  designed.  Therefore,  partial  sliding  mode 
control  results  in  the  stability  of  Xi  for  every  x2,  while  the 
conventional  technique  results  in  the  stability  of  all  states,  i.e.,  x. 

Consider  the  following  nonlinear  control  system  [15]: 

x=/(x,u),  x(t0)  =  x0  (5) 


and  greater  than  a  nonzero  lower  bound.  Under  this  assumption, 
the  partial  sliding  manifold  is  selected  as  s=£-(p{rj)=0,  where 
cp{rj)  is  a  continuously  differentiable  function  with  cp{ 0)=0.  In  this 
case,  the  sliding  manifold  is  a  function  of  only  the  first  part  of  states, 
i.e.,  £  and  rj.  It  was  shown  in  [15]  that  if  (p{rj )  designed  such  that: 

^^-fat(>lMil)<X2)<-y(\\ri\\)  V(f/,x2)  eRn'~k  x  R"2  (10) 

where  V  :  Rn'~k^>R  is  a  positive  definite,  continuously  differentiable 
function,  and  y(.)  is  a  class  K  function,  then  the  reduced-order 
model  of  Xi-subsystem,  i.e.,  rj  =fa\  (rj.cpirj)^)  is  asymptotically 
stable  with  respect  to  r\  for  all  x2  e  R"2  and  by  using  the  feedback 
control  law  (11),  only  trajectories  of  a  part  of  state  variables  (i.e.,  r\ 
and  0  are  forced  to  reach  the  partial  sliding  manifold  (s=£- 
(p(r])=0)  and  stay  on  it. 

u  =  (GkEkr'(^fa,-fb^J  (11) 


where  xeRn  is  the  state  vector  and  ueRm  is  the  input  vector.  Let 
vectors  Xi  and  x2  denote  the  partitions  of  the  state  vector, 
respectively.  Therefore,  x=  (x\,xl)T  where  X\  eRn\  x2eRni  and 
n\+n2  =  n.  As  a  result,  the  nonlinear  system  (5)  can  be  divided 
into  two  parts: 


Now,  consider  the  system  (9)  in  the  presence  of  uncertainties 
and  disturbances: 

n  =fa  1  (tlZX2) 

t  =/6l(f/.^X2)+gM(Xi,X2)U  +  /li(X,,X2,t)  (12) 


Xi (t)  =  F,  (Xi (t),x2(t),u),  X, (to)  =  x,o 

x2(t)  =  F2(Xi(t),x2(t),u),  x2(t0)  =  x20  (6) 

Definition  1.  The  nonlinear  control  system  (6)  is  said  to  be 
asymptotically  stabilizable  with  respect  to  x1  for  every  x2  e  Rn 2 ,  if 
there  exists  some  admissible  feedback  control  law  u=k(x i,x2),  which 
makes  the  nonlinear  system  (6)  asymptotically  stable  with  respect  to 
Xi.  It  means  that  in  the  closed-loop  system  for  every  s  >  0  and 
x20  e  Rni ,  there  exists  <5  =  <5(x2O)>0  such  that  ||xi0||<<5  implies 
|  |xi(t)|  |  <  s  for  all  t>  0  and  also  limt^ooXiO:)  =  0. 

Now,  suppose  the  k\ -subsystem  in  Eq.  (6)  is  affine  with  respect 
to  input  (the  x2 -subsystem  may  have  the  general  dynamical 
form).  Therefore, 

m 

Xi  =/1(Xi,X2)+  ^g,i(Xl,X2)Ui 
1  =  1 

x2(t)  =  F2(x1,x2,u)  (7) 

where  ut  is  the  ith  component  of  input  vector  u.  Also,  gu  e  Rn\  for 
i=l,2,...,m.  Let  k  be  the  number  ofgi*  vectors  which  are  nonzero. 
Flence,  k  is  the  number  of  components  of  input  vector  which 
appear  in  X\ -subsystem  and  thus  0  <  k  <  m.  Now,  we  assume  that 
k  #=  0.  By  augmenting  the  k  nonzero  vectors  gu  in  a  matrix,  the 
nonlinear  system  (7)  can  be  rewritten  as  follows: 

Xi  =/i(x1,x2)  +  G1(x,.x2)u 

x2  =  F2(x,,x2,u)  (8) 

where  UeRk  is  the  reduced  version  of  input  vector  u,  that  contains  k 
control  variables  appearing  in  X\ -subsystem,  Gi(xi,x2)  is  a  rq  x  k 
matrix  where  its  columns  are  the  k  nonzero  vectors  glim  In  this  case, 
the  task  is  to  find  an  appropriate  u,  which  guarantees  partial 
stabilization  of  nonlinear  system  (8)  with  respect  to  Xi  by  PSM 
technique. 

Let  Z:D-»Rni  be  a  diffeomorphism  that  transforms  only 
X\ -subsystem  in  (8)  into  the  regular  form.  Thus, 

'  (  d=fa\(n^2) 

<  =fbiQb&x2)+gbl(XuX2)u  (9) 

x2(t)  =  F2(xlfx2,u) 

where  [rf  £T]T=Z{x i),  £,ueRk  and  rjeRni~k.  Suppose  gbi(xi,x2)= 
G/c(xi,x2)£/c(xi,x2),  where  E/<(xi,x2)eR/<xk  is  a  nonsingular  matrix  and 
G/c(xi,x2)ef?/cx  k  is  a  diagonal  matrix  whose  elements  are  positive 


*2(0  =  F2  (Xi  ,x2  ,  u)  +  h2  (Xi  ,x2 , 0 


where  and  h2  are  unknown  nonlinear  vector  functions  which 
satisfy  the  matching  conditions;  that  is,  the  uncertain  terms  enter 
the  state  equations  at  the  same  point  as  the  control  input.  These 
terms  could  result  from  modeling  error  and/or  uncertainties  and 
disturbances,  which  exists  in  any  practical  problem.  In  a  typical 
situation,  hi  and  h2  are  unknown,  but  some  information  about 
them,  like  an  upper  bound  on  their  norm  may  be  known.  Suppose 
hi(xi,x2,t)  =  Bfc(Xi,x2,t)co(Xi,x2,t),  where  BkeRk*k  is  a  diagonal 
matrix  and  cdeRkxl  may  be  an  unknown  vector  with  a  known 
upper  bound  y(xi,x2,t),  i.e., 

|co/(x1,x2,t)|  <7,(xi,x2,t)  for  i  =  !,...,/<  (13) 


In  this  case,  the  task  is  to  design  u  to  force  s  toward  the 
manifold  s=0.  It  was  shown  in  [15]  that  if  ratio  BkiIGki  satisfies  the 
below  inequality  for  i=l,...,k, 


Bki(xux2,t) 

Gid(x  i,x2) 


<Pi(x  i,X2) 


(14) 


where  p{x ltx2)  is  a  known  positive  continuous  function  and  Bki 
and  Gki  are  the  ith  elements  of  diagonal  matrixes  Bk  and  Gk  on  its 
diagonal.  Then  the  following  feedback  law  guarantees  robust 
partial  stabilization: 

u=£r1{cr,(^/ai-/M)+^}  as) 

where  v  e  Rk  and  its  ith  component  is 

Vi  =  -ft(xlfx2)sgn(Si)  (16) 

and  ft(xltx2)  > pf(x1,x2)yi(x1,x2,t)  +  jS0i  where  poi  is  a  positive 
constant.  Using  (15)  guarantees  that  all  trajectories  starting  off 
the  manifold  s=0  reach  it  in  a  finite  time  and  those  on  the 
manifold  cannot  leave  it.  Restricting  the  motion  on  this  manifold 
guarantees  stability  with  respect  to  £  and  r\  (partial  stability)  for 
allx2  eR"2  [15]. 


Remark  1.  For  the  case  0  <  k  <  m,  there  are  l=m-  k  components  of 
control  input  vector  u  not  appearing  in  the  x,\  -subsystem  and  there¬ 
fore  do  not  appear  in  u,  which  appear  in  the  x2 -subsystem.  By  putting 
these  components  in  vector  u  e  Rl,  this  vector  may  be  chosen  in  a  way 
that  the  response  of  x2  satisfies  the  constraints  related  to  practical 
considerations  of  the  application.  If  there  is  an  uncertain  term 
h2(.X]jc2,t')  in  the  x2 -subsystem,  the  additional  term  v  could  be 
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designed  in  a  way  that  the  control  law  u  +  v,  guarantees  the  desirable 
behavior  for  x2  in  the  presence  of  h2(x ltx2,t)  (as  illustrated  in  the  next 
section). 

Remark  2.  Since  discontinuous  controllers  suffer  from  chattering, 
one  way  to  alleviate  this  problem  is  to  consider  an  approximation 
of  the  signum  function  by  a  saturation  function  with  a  high  slope. 

4.  Guidance  law  design 

In  this  section,  we  employ  the  PSM  technique  to  fulfill  the 
design  task  of  the  guidance  law. 

4.1.  PSM  guidance  law 


which  implies  exponential  decaying  of  rj.  Furthermore,  according  to 
Eq.  (15): 

v=E*1{(jinfal-fbl)+v} 

=  — G^(— M£-/m)—  v  (21) 

where  Ek  is  substituted  from  (18b).  A1so,/m  is  given  in  (18a).  The  final 
step  is  to  construct  an  extra  control  v  that  can  compensate  for  the 
effect  of  target  acceleration  vector  (external  disturbance).  For  this 
purpose,  assume  that  the  components  of  acceleration  vector  of  the 
target  satisfy  the  following  assumption: 

Assumption  1.  There  exists  a  vector  function. 

.*2)  =  B»i  ,y2]T  =  [ye(*i  ,x2),y^(xi  ,x2)]r  (22) 


Let  x-t  =  [t]T  <ff,  where  =  f/2]T=[0  4>V  and  £  =  £2]T  = 

[9  4>]T.  Also,  x2  =  [x2i  x22]  =  [r  r]T.  co=\cor,me,co^)T  and  u={ur,UeM^)T 
(the  acceleration  vectors  of  target  and  missile),  are  assumed  as 
external  disturbance  and  input  control  vector,  respectively.  The 
kinematics  equations  (2)  can  be  rewritten  in  the  following  nonlinear 
state-space  equation: 


Xi  -subsystem 


r«7i  =  {i 
'72  =  52 

<  ^=2^2tanri2-2^—^njiue+^^(oe 
j2  =  -2  ^  -  d  cos  t]2  sin  t]2  -  i  +  J- 


x2 -subsystem 


J  X21  =  x22 

f  X22  =  X21  £2  +X21  5i  COS2  t]2-Ur  +  CO r 


(17) 


Note  that,  the  X\ -subsystem  in  dynamical  system  (17)  is  already 
in  the  regular  form.  By  substituting  (17)  in  the  structure  (12),  one 
may  define: 


fa\  ~ 


& 


fb\ 


2{1{2tanff2-2**iC 
cosr]2  sinil2 


(18a) 


1  0  1 

*21  COS  7/ 2 

-1 

0  ■ 

r  1  0 1 

X21  cos  77  2 

;  u  = 

~Uq~ 

1 - 

O 

1 

?|- 

1 _ 

_  0 

-1 

0  -L 

*21 

-e- 

3 

(18b) 


hi(x1,x2,t)  = 


X21  cost; 2  ^0 


r  1 

0  " 

~coe~ 

*21  cos  77  2 

0 

*21 

CO 0 

Bk 

CO 

(18c) 


and 


x22 

'  0  ' 

O' 

X2lf2+X2lf?  COS2  t]2~Ur 

,  h2(x1,x2,t)  = 

COr 

— 

1 

C Or 

08. d) 


It  is  worth  noting  that  throughout  the  interception  process,  it  is 
logical  to  assume  that  x2i  —  r  >  0  and  1 </> 1  <  7c/2  [5].  Also,  x2i  is  smaller 
than  its  initial  value,  i.e.,  x2i(0)=ro.  This  implies  that  Gn  >  l/r0  >  0  for 
z=l,2,....  To  fulfill  the  PSM  technique,  we  choose  the  partial  sliding 
surface  as 

S  =  (SUS2)T  =  Z+Mri  (19) 

where  cp{r])=-Mr]  and  MeR 2x2  is  a  positive  definite  matrix.  It  is 
worth  nothing  that  if  i-rj  lies  on  the  partial  sliding  surface,  the 
reduced  model  have  the  form: 


(20) 


where  ye  and  y $  are  non-negative  functions  such  that  |cof  |  <  y,  for 
i=l,2. 

Furthermore,  \BkiIGki\  =  \  for  i = 1 ,2  (refer  to  (18b)  and  (18c)). 
Therefore,  p/(xi,x2)  =  l.  Now,  by  considering, 

vi  =  -A(x1,x2)sgn  (s^  (23) 

where 

ft(x  1  ,x2)  >  pi(x  1  ,x2)7j(xi  ,x2)  +  poi 

=  y,(xi  ,x2) + Poi  (24) 

As  a  result,  the  reduced  control  vector  (21)  achieves  the 
desired  performance  for  Xi.  Now,  in  order  to  complete  the  design 
process,  u  =  ur  should  be  obtained  to  attain  the  appropriate 
behavior  of  x2  (refer  to  Remark  1).  Assume  |cor|  <yr(xi,*2)  and 
take 

ur=x2\il+x2\i\  cos2  t]2-ax22  +  Vr  (25) 

where  a  >  0  and  i/r=yr(x i,x2).  By  substituting  (25)  in  (17),  one  has 
X22=(JX22-yr  +  (Dr 
Thus, 


x22(t)=x22(0)erTt+  [  (— yr+cur)e0'(t_T)  dr 
Jo 

<x22(0)eat-  [  yreCT(t_T)  dx  +  f  yrea(t~T)  dz  (26) 

Jo  Jo 

Consequently,  x22(t)  <x22(0)  <  0.  Therefore,  by  choosing 
-C=x22(0),  Proposition  1  is  satisfied  and  the  desirable  behavior 
for  x2  is  achieved.  In  addition,  the  relative  distance  will  turn  into 
zero  within  a  finite  time.  Clearly,  higher  values  of  a  makes  a 
shorter  time  of  interception;  however,  its  adjustment  should  be 
done  with  respect  to  physical  limitations  [14]. 

Consequently,  the  PSM  guidance  law  (27)  guarantees  inter¬ 
cepting  of  maneuvering  targets  within  a  finite  interception  and 
zero-miss  distance. 


<5’ (-«{-/«)-■>  (27) 

[  ur  =  x2 i£2+X2i5?  cos2  rj2-ax22+yr 

where  £  =  [£,  £2]T  =  [()  <£]T,  r}  =  [r\x  r]2]T=[9  4 )]T .  s=(sif2)T=£+ 
Mr],  (MeR2  x 2  is  a  positive  definite  matrix),  x2  =  [x21  x22]  =  [r  r]r, 
and  v  =  [Ti,I72]rwhere  =  -^(Xi,x2)sgn(Sj),for  i=l,2.  Also,  Gk  and 
fbx  are  given  in  Eqs.  (18a)  and  (18b),  respectively. 


4.2.  Computer  simulations 

Numerical  simulations  are  performed  to  illustrate  the  effec¬ 
tiveness  of  the  PSM  guidance  law  against  a  highly  maneuvering 


p  +Mr]  =  0 
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target.  The  target  is  assumed  maneuvering  at  the  following 
trajectory: 

co(t)  =  70  sin(0.5t)7fr  +  70  sin(O.5t+7r/4)"?0  +  7O  cos(0.5 t)~e  $ 

(28) 

The  initial  states  are  chosen  as  r0=5km,  0o=7t/3,  0o=7t/3, 
0  =  0.08  rad/s,  </>  =0.06  rad/s  and  r=  -300  m/s.  The  initial  condi¬ 
tions  in  the  Inertia-frame  of  the  missile  and  the  target  are  given  by 

•  For  the  missile:  xM(0)=yM(0)  =  0  m 

•  For  the  target:  xr(0)  =  5000  cos 0O  cos 90m,  yT( 0)  = 

5000  cos0o  sm60m,  zT( 0)  =  5000  sin 4>0m  and 

vTx(0)=vTy(0)=VTz{0)=W0  m/s. 


In  addition,  /?Oi=/?o2  =  0.1,  cr  =  0.03  and  the  three  non-negative 
functions  yr(xi,x2),  ye{x i,x2)  and  y^(x i,x2)  are  assumed  as 
yr=y0=')4  =  7O.  Also,  M=diag{2.41,  2.41}  is  selected.  In  the  mean¬ 
time,  for  the  purpose  of  attention  of  chattering,  instead  of  sgn(s/), 
the  saturation  function  sat(s//e)  with  e=l,  is  used. 

The  PSM  guidance  law  is  compared  with  the  SM  guidance  law 
presented  in  [5].  It  is  because  that  the  same  state  space  model  has 
been  used  in  [5].  Additionally  SM  guidance  law  has  been  compared 
with  some  existing  guidance  laws  in  [5]  and  has  achieved  better 
results.  It  is  shown  in  this  section  that  the  proposed  guidance  law 
has  gained  the  better  performance  in  terms  of  interception  time  and 
control  effort,  compared  to  the  SM  guidance  law.  In  SM  guidance 
law  stability  of  all  states  is  under  consideration  in  contrast  to  the 
PSM  guidance  law.  Numerical  simulations  are  given  in  Figs.  2-4. 
As  seen,  the  PSM  law  has  totally  gained  the  better  performance  in 
terms  of  interception  time  (tfi  and  control  efforts. 

Fig.  2  shows  the  trajectories  of  the  missile  (M)  and  the  target 
(T)  in  the  Inertia-frame,  where,  Cl  and  C2  are  the  collision  points 
for  the  PSM  guidance  law  and  the  SM  guidance  law,  respectively. 
The  interception  time  is  9.04  s  for  the  PSM  guidance  law  and 
12.38  s  for  the  SM  guidance  law. 

Fig.  3  illustrates  the  guidance  commands.  Fig.  3(b)  and 
(c)  show  that  the  PSM  guidance  law  requires  fewer  efforts  in 
tangential  components  of  input  vector,  i.e.,  ue  and  14,  than  those 
for  the  SM  guidance  law.  Also,  Fig.  3(a)  demonstrates  that  the 
amount  of  changes  of  radial  control  effort  for  the  PSM  law  is  less 
than  that  of  SM  guidance  law.  In  brief,  the  control  energy 
(Jq7  uTudt)  of  the  PSM  law  is  much  ur  less  than  that  of  the  SM  law. 


C 

2500 


- Sliding  Mode  Law 

- Partial  Sliding  Mode  Law 


5.  Conclusion 

In  this  paper  a  new  viewpoint  to  the  guidance  problem  was 
presented.  It  was  demonstrated  that  in  practical  approach  to 


_500  - 1 - 1 - 1 - 1 - ' - ' - 1 - 1 - 1 - 1 - ' - 1 - 

01  2  34  56  7  8  9  10  11  12  13 

time  ( s ) 

Fig.  3.  (a)  Radial  guidance  command  (ur),  (b)  tangential  guidance  command  tie 
and  (c)  tangential  guidance  command  u ^ 

guidance  problem,  the  desirable  behaviors  of  state  variables  are 
different  with  respect  to  each  other,  and  the  asymptotic  conver¬ 
gence  behavior  is  not  ideal  for  all  state  variables.  Therefore,  based 
on  PSM  technique,  a  new  robust  nonlinear  guidance  law  was 
developed  and  its  effectiveness  in  interception  of  maneuvering 
targets  was  shown  analytically  and  also  through  simulations.  In 
the  simulation  section,  complicated  and  highly  maneuvering 
targets  (sinusoidal  targets)  was  considered  and  it  was  shown  that 
the  proposed  method,  in  comparison  to  sliding  mode  method,  can 
lead  to  better  results.  It  is  because  this  method  (in  contrast  to 
other  robust  nonlinear  guidance  laws  like  SM  guidance  law)  does 
not  need  to  regulate  r  to  a  negative  constant  c.  The  value  of  this 
constant,  considerably  affects  the  performance  of  guidance  law 
and  in  order  to  determine  it  efficiently,  the  maneuvering  model  of 
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a  b 


Fig.  4.  (a)  Relative  distances  between  the  missile  and  the  target  (r),  (b)  radial  component  of  relative  speed  (r),  (c)  pitch  LOS-angular  rate  (6)  and  (d)  yaw  LOS-angular  rate  ($). 


the  target  should  be  known  a  priori  for  the  guidance  system, 
while  it  is  not  always  possible  in  practice. 
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